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Abstract: The Hubble constant during the last stages of inflation in a broad class of mod- 
els based on the KKLT mechanism should be smaller than the gravitino mass, H < m^ii- 
We point out that in the models with large volume of compactification the corresponding 
constraint typically is even stronger, H < m ^/2i 111 Planck units. In order to address this 
problem, we propose a class of models with large volume of compactification where infla- 
tion may occur exponentially far away from the present vacuum state. In these models, 
the Hubble constant during inflation can be many orders of magnitude greater than the 
gravitino mass. We introduce a toy model describing this scenario, and discuss its strengths 
and weaknesses. 
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1. Introduction: Inflation and the Gravitino Mass 

A realistic model of the universe should be consistent with all constraints, not just from high 
energy physics experiments but also from cosmological observations. Inflation is the leading 
scenario for the early universe and low-energy supersymmetry is the leading proposal for 
new physics subject to experimental test at the LHC. It is therefore important to develop 
models capable of including both these desirable features of physics beyond the standard 
model. 

Many recent attempts to implement inflation in the context of string theory are based 
on the KKLT mechanism of vacuum stabilization [1] and its generalizations [2]. Until very 
recently, the standard assumption of many of these inflationary models was that inflation is 
a very high energy scale phenomenon, and therefore one can construct inflationary models 
quite independently of the requirements of the low-energy SUSY phenomenology. 

However, recent studies of inflationary models in string theory revealed a rather unex- 
pected fact: in the simplest models based on the KKLT mechanism the Hubble constant 
H should be smaller than the present value of the gravitino mass [3], 



The reason for this bound is that the gravitino mass at the supersymmetric KKLT min- 
imum, with DW = before the uplifting, is given by 3m5y 2 = (V^dsl- Uplifting of the 
AdS minimum to the present nearly Minkowski vacuum occurs by adding to the potential 
a term of the form C/a n , where a is the volume modulus and n = 3 for a generic com- 
pactification and n = 2 for the highly warped throat geometry. Since the uplifting is less 
significant at large a, the energy barrier to decompactification created by the uplifting is 
generically slightly smaller than |Va<is|: Ub aiT ier < IK^S 1 ! ~ 3m?, 2 . However, the volume 
modulus couples to all sources of energy due to the Weyl rescaling always present when 
deriving the 4-dimensional action. In particular, the energy of the inflaton field 4> will give 
an additional uplifting of a similar type: AV((p, a) ~ V((p)/cr n . As this is also proportional 
to an inverse power of the volume modulus, it is larger at the minimum of the KKLT 
potential than at the top of the barrier. Therefore adding a large vacuum energy density 
to the KKLT potential, as required for inflation, may destabilize the minimum by uplifting 



H < m- 



'3/2 • 



(1.1) 
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it to a height greater than the height of the barrier, see Fig. [l| In typical KKLT-type 
models this leads to vacuum destabilization if the added energy density V((f>)/a n , which is 
responsible for inflation, is much greater than the height of the barrier V\> a 
Since H 2 



< "\m 2 M 2 
~ om 3/2 iV1 P- 



AV(</>, cr)/3, this leads to the bound (LI) (see [3] for a more detailed discussion 
of this issue, while a similar problem in a slightly different context was also found in [4]). 
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Figure 1: The lowest curve with dS minimum is the potential of the KKLT model. The second 
one shows what happens to the volume modulus potential when the inflaton potential Vi n fl = 
added to the KKLT potential. The top curve shows that when the inflaton potential becomes too 
large, the barrier disappears, and the internal space decompactifies. This explains the constraint 
H < to 3 / 2 . 



In KKLT-based models, it therefore seems that for a gravitino mass ~ ITeV the 
Hubble constant during the last stages of a string theory inflation model should be quite 
low, H < 1 TeV, which is ten orders of magnitude below the often discussed GUT inflation 
scale. Therefore if one believes in standard SUSY phenomenology with m 3 / 2 ^ 0(1) TeV, 
one should find a realistic particle physics model where the nonperturbative string theory 
dynamics occurs at the LHC scale or even lower (the mass of the volume modulus in the 
KKLT scenario typically is not much greater than the gravitino mass) , and inflation occurs 
at a density at least 30 orders of magnitude below the Planck energy density [3]. For a 
recent analysis of this issue see e.g. [5] and for a discussion in the context of the heterotic 
string see [6]. 

This problem is quite generic. For example, recently a new interesting mechanism of 
moduli stabilization was proposed, which is based on the models with compacification on Nil 
manifolds with negative curvature [7]. This mechanism presents a significant modification 
of the compactifications on flat Calabi-Yau spaces, as suggested by the assumption of the 
low scale supersymmetry. And yet, the same constraint H < w-3/2 remains valid for the 
inflationary models in this scenario [8]. 

The situation becomes even trickier in the large volume models of vacuum stabilization 
[2]. In such models the height of the barrier is much smaller, Vbamer ~ m \/2^ p - ^ n ^ ms 
case, the constraint that the inflaton potential should not be much greater than the height 
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of the barrier leads to the bound (in units M p = 1) 

H < m^l . (1.2) 

For 771.3/2 ~ 1 TeV this inequality implies that the Hubble constant during inflation in this 
class of models [9, 10] cannot exceed O(l) KeV, which is an extremely strong constraint. 

There do exist proposals of low-scale inflationary models, for example the so-called 
MSSM inflation, which may occur for H ~ 10 GeV or even for H ~ 10 MeV [11]. Ref. [12] 
also contains a discussion of models where inflation may occur at extremely low scales, with 
an example of a model for which H ~ 10~ 7 eV. In particular, if the inflaton potential energy 
at H ~ 1 KeV could instantly transfer to thermal energy, the corresponding temperature 
would be about 10 6 GeV, which is much greater than the critical temperature of the 
phase transition in the standard model. If this instantaneous transition is achievable, the 
temperatures would then be sufficiently high for the subsequent generation of a baryon 
asymmetry. 

One can find models with a very low-scale inflation in the context of the KKLT or 
large volume scenarios, since the energy scale is exponentially sensitive to the parameter 
a of the nonperturbative superpotential W = Wq + Ae~ aT [1]. However, models of this 
type are very non-traditional, and their parameters are substantially different from the 
parameters of all current existing models of string theory inflation. Furthermore, as the 
required value of the slow-roll epsilon parameter is given by e ~ (£ inf /6 x 10 16 GeV) 4 , low- 
scale inflation substantially increases the amount of fine-tuning required in the inflaton 
potential. It is important to know whether this tension between high-scale inflation and 
TeV supersymmetry is unavoidable or whether it is simply a consequence of the assumptions 
used so far in inflationary model-building. 

This is not the first time that string theory and supergravity have encountered cosmo- 
logical problems associated with the small value of the gravitino mass and of the moduli 
fields. The famous gravitino problem and the cosmological moduli problem are haunting 
us for more than two decades [13-15]. Now we see that the smallness of the gravitino 
mass leads to an additional problem in the context of string cosmology [3,4]. This problem 
would disappear if one would consider supersymmetric models with large gravitino mass, 
for example [16, 17], or used a solution to the hierarchy problem different to that of TeV 
supersymmetry. 1 

There exist ways to address this problem without increasing the value of the gravitino 
mass. For example, one may consider KKLT models with the racetrack superpotential 
containing at least two exponents and find parameters such that the supersymmetric min- 
imum of the potential, even prior to uplifting, occurs at zero energy density [3], which 
would mean 1113/2 = 0. By a slight change of parameters in this class of models, which are 
sometimes called KL models, one can get a gravitino mass that is nonzero but still much 
smaller than the height of the barrier, removing the constraint H < m 3 / 2 . In particular, 

1 For other problems with high values of the Hubble constant in string inflation see [18]. 
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one can use the KKLMMT brane inflation model [19-25] and implement it in the context 
of the KL scenario with H 3> m 3 / 2 . 

The difficulty with this solution is that if we want to have H many orders of magnitude 
greater than m 3 / 2 ~ 1 TeV, we need to fine-tune the parameters of the model to a corre- 
sponding accuracy. The origin of the electroweak scale is then a kind of accident, reducing 
the attractiveness of super symmetric solutions to the hierarchy problem. However, this 
class of models has certain advantages from the point of view of vacuum stabilization [3], 
so it might happen that the required fine-tuning is not unreasonable. 

Another possible solution of this problem was recently proposed in [26] in the context 
of the volume modulus inflation in the KKLT scenario. This model extended the KL 
model to involve triple gaugino condensation in the superpotential and also required a 
modification of the Kahler potential. It introduces 6 new parameters (two real parameters 
and two complex ones), and fine-tuned three of these parameters with accuracy ranging 
from 10~ 4 to 1CT 7 . This clearly demonstrates that it is quite difficult to avoid the constraint 
H < m 3 /2 in such models, but nevertheless it is encouraging that it is possible to do so. 

In this paper we will concentrate on the models with large volume of compactification 
[2] and propose another possible resolution of the gravitino mass problem, aiming at making 
a transition of scales from E^ n j to E susy natural. We shall discuss our idea and certain 
issues associated with it in the next sections. 

2. Disentangling H and m 3 / 2 : the basic idea 

The idea we revisit in this paper is the assumption used so far in building supergravity 
models of inflation that the true minimum of the scalar potential is relatively close in field 
space to the locus at which inflation ends. We instead propose that inflation should end 
with a runaway in field space, with the true minimum lying a very long distance in field 
space from the location where inflation occurs (with the model discussed below, the distance 
in canonically normalized field space will correspond to approximately twenty Planckian 
distances) . The problem we are trying to address is that characteristic inflationary energy 
scales are much larger than those appropriate for supersymmetry breaking. The advantage 
of a runaway epoch is that evolution along runaway potentials (e.g. of the form V((p) ~ 
y^ e -\<f>^ - g Qne Q £ £ ew e ffi c i en t ways of naturally dissipating large quantities of energy 
and reducing the scale of the potential by many orders of magnitude. If the true minimum 
of the scalar potential lies a long way along the runaway direction, it naturally has much 
lower characteristic energy scales than apply during inflation. 

In this case supergravity models of high-scale inflation consistent with low-scale super- 
symmetry breaking, m 3 / 2 < 1 TeV in vacuo, should have three stages. In the first, inflation 
occurs at high energy scales with m 3 / 2 S> 1 TeV during inflation. In the second, inflation 
ends with the fields fast-rolling towards a runaway direction. For example, in the model 
below this is due to inflation occurring near an inflection point in the volume direction. 
As inflation ends it is necessary that trace quantities of radiation be generated to act as a 
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seed for an attractor solution. In the third stage, the presence of small initial quantities 
of radiation drives the fields to an attractor solution. The attractor solution applies dur- 
ing the runaway epoch and dissipates energy. The scaling nature of the attractor solution 
avoids overshooting and guides the fields into the global minimum of the potential in which 
m 3/2 ~ 1 TeV. 2 This scenario is illustrated in Fig. g. 

The justification for the existence of a minimum at very large values of the volume, 
far along the runaway direction, is the large volume scenario [2], where the inclusion of o! 
corrections into the KKLT framework generates a new minimum of the scalar potential at 
exponentially large values of the volume, with hierarchically small values of 7713/2- 

V 




Figure 2: An illustration of the scenario put forward in this article. At relatively small vol- 
ume, high-scale inflation occurs due to fine-tuned quantum corrections. After inflation the volume 
modulus evolves over a long range of many Planck scales, eventually settling in the large volume 
minimum with TeV gravitino mass. Although the barrier protecting from decompactification is 
very small compared to the initial energies, an attractor solution guides the fields to the minimum 
and prevents overshooting. 

To illustrate this idea, we start by studying moduli evolution in the following toy model 
describing a field with a potential 

V = y ((l-e$ 3/2 )e~^ # +Ce -io*Ms + £ e -ii*/V6 + Se -VE^ (21) 

The particular form of this potential is motivated by that arising as the effective potential 
for the volume modulus in the large volume models. The connection to the large volume 

2 For a recent discussion in the context of M-theory compactifications of the overshooting problem and 
how to avoid it, see [27]. 
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models and the supergravity origin of the above potential will be discussed in the next 



section. The 1-modulus potential (2.1) will not represent a complete model but will allow 
us to capture several key features of our proposal. 

In (|2.ip $ is the canonically normalized volume modulus, = y/3/2 log t\, with V = 



3/2 

r b , so Tb is the volume of a 4-cycle. The first two terms of the potential correspond to 
the effective F-term potential for the volume modulus in the large volume models. The 
structure of these terms generates a minimum at large values of ~ e -2 / 3 . The 

definition of $ as log V implies that this is equivalent to the existence of a minimum 
at exponentially large volumes. As in string theory the gravitino mass is given by = 
Mp(Wo/V), where Wo is the constant (flux) superpotential, hierarchically large volumes 
corresponds to hierarchically low gravitino mass. The final term, depending on 5, is the 
uplift term needed to lift the minimum of the potential to Minkowski space. 

The second and third terms, depending on C and D, are only important at small 
volumes. These terms are not derived from string theory but are included in a phenomeno- 
logical fashion in order to generate an inflection point, and thus inflation, at small volume. 
Besides the large volume minimum, occurring for large values of <J> for which the C and 
D terms are irrelevant, by tuning C and D 3 this sum of exponentials can have a non- 
monotonic shape for relatively small values of the field <£. 

We will consider the following set of parameters, which are determined by the require- 
ment of obtaining inflation satisfying the constraint which follow from the existing obser- 
vational data: V = 1.45 x 1(T 14 , e = 0.013, C = -3, D = 2.3045, 5 = 0.06155 x 10~ 10 . 
The shape of the potential for this parameters at small 3? is shown in Fig. [|. As we see, 
it has an inflection point at ~ 1.3, where inflation takes place. The Hubble constant at 
that time is about 10~ 9 M P ~ 10 9 GeV. 

The potentials of this type have been used by many authors, see e.g. [11,21-25,28-30]. 
The typical feature of this class of models is that by tuning their parameters, the potential 
can always be represented as a sum of a linear term and a cubic term in a vicinity of the 
inflection point, 

V = V($ ) (l + Ai(3> - $ ) + y (* - $o) 3 + • (2.2) 

Here 3>o is the position of the inflection point. (In our case, $0 ~ 1-3, in Planck units.) The 
numerical values of Aj are determined by fine-tuning of the values of the parameters in our 
full expression for V. By making the linear term as small as possible, one can maximize 
the exponential growth of the universe during inflation, which may serve as a possible 
justification of the fine-tuning of the potential required for a long stage of inflation [29]. 

In the case where inflation is dominated by the cubic term, inflation starting at the 
inflection point is eternal, the spectral index n s ~ 0.93, and the amplitude of perturbations 



3 Notice that in full string compactifications, corrections to the Kahler potential such as may generate the 
parameters C, D will be functions of both fluxes and complex structure moduli and so would be expected 
to take on many possible values in the landscape. 
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of metric produced during inflation satisfies the COBE-WMAP normalization for 



(2.3) 



in units of Planck density [29]. We calculated V(3>o) an d M an d tuned the parameters 
of our potential to satisfy this constraint and make Ai vanishingly small. (If one reduces 
the degree of fine-tuning and considers a theory with a non- vanishing negative Ai, one can 
significantly increase n s , but it will simultaneously decrease the degree of expansion of the 
universe during inflation [29].) 

V 
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Figure 3: The potential at small <&. Inflation occurs near the inflection point at $ 
units. 



1.3, in Planck 



After inflation, the field enters a runaway period until it is captured by the minimum 
at <3? ~ 19, see Fig. ||. The parameters of our model were tuned to make the vacuum 
energy nearly zero in the minimum, and to have the gravitino mass there in the TeV range, 
m 3/2 = 0(l)TeV. This last statement assumes that the gravitino mass is determined 
entirely by the volume, taking Wo ~ 1 as usually done in the large volume models of 
moduli stabilisation. 

V 
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Figure 4: The potential at large <&. Vacuum state corresponds to the minimum at $ ~ 19. 



Let us study the cosmological dynamics of rolling moduli in the potential (2.1). The 
C and D terms are only important at small $ and are highly suppressed for large <!>. 
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Conversely, the 5 and e terms only become important at large <£. These generate the 
minimum at exponentially large volumes and uplift it to Minkowski space. However, for 
smaller values of the volume these terms contribute highly subleading corrections to the 
potential. The upshot is that for a very large range of the volume the dominant term in 



the potential (2.1) is simply 



V = V e V 2 * (2.4) 

with <I> canonically normalized. In this regime the cosmological dynamics reduces to that 
of a pure exponential potential. 

Evolution along runaway directions of this type is associated with the overshoot prob- 
lem [31] and it is necessary that the fields are able to locate the global minimum of the 
potential. In the presence of any additional sources (such as radiation) exponential poten- 
tials have attractor solutions in which the different components of energy track the total 
energy as a constant fraction. This attractor solution can guide the fields to the global 
minimum without overshooting. In order to avoid decompactification and the overshooting 
problem, it is necessary that a small component of radiation be present in the evolution 
equations. We study the evolution with radiation present immediately at the end of infla- 
tion. 4 



The evolution of fields in pure exponential potentials (2.4) has been studied for a long 
time. Naively <3? will rapidly roll down the exponential slope into a kination phase. However 
in the presence of a background fluid, such as radiation, this does not occur [32-39]. Since 
kinetic energy density falls faster than density of radiation, 



1 1 

Prad ~ ^J, PRE ~ ^ 



the radiation comes to dominate the energy density, contributing additional Hubble friction 
to the field dynamics. The late-time attractor is a scaling solution in which the different 
components (radiation, potential energy, kinetic energy) constitute fixed fractions of the 
overall energy density. 

Following [32,34], the equations of motion are conveniently formulated in terms of the 
variables x and y defined by 




y 



x 2 and y 2 equal tl$ t kin and 0$ )PO £ respectively. It is also convenient to use N = In a as the 
time variable. The equations of motion are 

x'(N) = -3s - I/ 2 + \x (2x 2 + 7 (1 -x 2 - y 2 )) , (2.5) 



4 A deficiency of the 1-modulus potential is that there is no way to generate this radiation, as any 
primordial radiation is diluted during inflation, and thus it has to be introduced by hand. With more 
complicated models, for example the 2-modulus case studied in the next section, there are ways to generate 
radiation at the end of inflation. 
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y'(N) = ^^y + \v (2x 2 + 7(1 -x 2 - y 2 )) , (2.6) 

H'(N) = -^(2x 2 + 7 (l-x 2 -y 2 )), (2-7) 

$'(JV) = V^x. (2.8) 

7 = 4/3 for radiation and 1 for matter. For an exponential potential V = Voexp(— A<3?) 
with A 2 > 6, there exists a scaling solution given by [32,33,35] 

2 _ 3 7 2 , _ 3(2 - 7 ) 7 3 7 

2A2' V ~ 2A 2 ' ~~ A 2 "' 

For the case at hand of A = yf^ and 7 = §, the attractor solution is 



^ 7= 2~7' = ^ = 81' ^ y=y2 = il' (2 ' 9) 



In the presence of any initial radiation or matter, this solution is the late-time attractor. 

Considering the full potential ( |2.1| ) , the attractor solution will exist so long as the pure 
exponential provides a good approximation to the potential. At small <&, this approxima- 
tion will hold soon after inflation has ended, once the C and D terms become negligible. 
Inflation occurs at V ~ 10 -18 , corresponding to $ ~ 1.3. At large the attractor solu- 
tion disappears as $ approaches the susy breaking vacuum and additional terms become 
important. We want the vacuum to satisfy m 3 / 2 ~ ITeV, and as the characteristic scale of 
the large-volume potential is is V ~ m^ 2 ~ 10~ 47 (assuming Wo ~ 0(1)), this is attained 
when $ ~ 19, see Fig. ||. 

The regime in which the potential is well described by ( |2.4| ) and the attractor solution 
is valid is therefore 

1-3 < * < 19, (2.10) 

a range of almost twenty Planckian distances. We note that despite the substantially 
trans-Planckian field range, the potential is under good control: the high-scale theory is 
understood and the potential is simply a decompactification potential. 

The physical minimum exists at $ ~ 19, shortly followed by a local maximum repre- 
senting the barrier to decompactification. In this regime the attractor solution is no longer 
present. It is necessary that the fields fall into the minimum rather than passing over the 
barrier and running away to decompactification. Whether this occurs or not depends on 
whether or not the field $ has located the attractor solution while in the range (2.10). If 



the attractor solution has been found, overshooting does not occur and $ settles into its 
minimum. This is illustrated in figure ||, which shows the behavior of the attractor solution 
as it approaches the minimum. 

Whether or not the fields locate the attractor solution depends on the initial conditions 
and in particular the initial amount of radiation present. This is illustrated in figure |6|. 
The initial conditions were chosen to be &i n it = 2 ~ ^ = with varying initial values 
of As expected, the larger the initial value of 7 the more rapidly the attractor 
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Figure 5: The evolution of the radiation-background attractor solution as it approaches the mini- 
mum. The solid dashed horizontal line shows the location of the barrier to decompactification, and 
the narrow horizontal line the location of the true minimum. The attractor solution settles at the 
minimum and does not overshoot. The different paths correspond to different initial conditions. 
N = In a is the time variable. 
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Figure 6: The evolution of the field with time for radiation (a) and matter (b) backgrounds. In 
both cases the field evolution is taken to start at &i n it = 2 with <&i n it = 0. In figure (a), the field 
evolution is shown for initial radiation densities f2 7 = 10~ 7 , 10~ 6 , 10~ 4 , 10~ 2 , 10 _1 . For £1 7 = 1CP 7 
or 10 -6 the field overshoots the minimum, whereas for the other values of 7 the field finds the 
tracking attractor solution and settles in the minimum. Figure (b) is for a matter background, with 
initial values of fl m = HT 11 , 1CT 9 , 1(T 7 , 1(T 5 , 1(T 3 . For n m = 10" 11 or 1(T 9 the field overshoots 
the minimum, whereas for all the other values of f2 m the field finds the tracking attractor solution 
and settles in the minimum. 



solution is found. As long as the attractor solution is found before the field reaches the 
decompactification barrier, overshooting does not occur. It is seen that even very small 
initial values of f2 7 are sufficient to avoid overshooting, with even O-y j n jt = 10 -4 being 
sufficient to prevent overshooting. This can be compared with the results of [34, 37] (see 
in particular figure 3 of [34]), where a similar question was studied in the context of a 
KKLT model, and to avoid overshooting J7 7 > 0.5 was required across the whole range of 
parameter space. 
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In the case studied here it is much easier to avoid overshooting, which can be achieved 
with Flyjnit C 1. It is easy to understand why this model is more efficient at avoiding 
overshooting than the models studied in [1,34] . If the runaway direction comes from gaugino 
condensation in either the dilaton or volume directions, then the canonically normalized 
potential is a double exponential e _e * (as the potential is ~ e~ T and the canonically 
normalized field is $ ~ InT). This is much steeper than the single exponential (e - *) 
present here, and so it is much harder to avoid overshooting. 

The above one-field model has illustrated the required form of the potential, with a 
region suitable for inflation at small volume, a minimum at very large values of the volume, 
and a region in-between where the potential is described by a runaway. In order to locate 
the attractor solution that will guide the field to the global minimum, it is necessary that 
a small amount of matter or radiation can be generated as runaway starts. An attractive 
feature is that the quantity of radiation required for the fields not to overshoot is very small 
($7 7 ~ 10~ 4 ). However, any primordial radiation will be diluted away during inflation. It 
is therefore necessary to extend the 1-modulus model to describe (partial) reheating and 
the generation of radiation at the end of inflation, as this is necessary to act as a source 
for the attractor solution. 



3. Two-Field Inflation 

We now consider a 2-field model which will justify the form of the potential used in the 
discussion of the 1-field model. The model we use is the following N = 1 supergravity 
theory 

A- = -21„(v + ? + T £j + ^), 

W = W + Ae~ asTs 

V = e K (k^D.WD-^W - 3|VF| 2 ) . (3.1) 

Except for the parameters C and D, this is the supergravity theory that describes the 
large- volume models for compactifications on the Calabi-Yau manifold PA 1 1 6 Q i [2]. if, 
and r s are Kahler moduli: t\j controls the size of the Calabi-Yau, whereas r s is the size of 
a small blow-up 'hole' in the Calabi-Yau. 

Although all input parameters are of order unity, this theory admits a vacuum at expo- 
nentially large values of the volume with hierarchically small supersymmetry breaking [2]. 
In fl3,ip additional corrections parametrized by C and D have been included. These are 
motivated by the existence of higher a' corrections to the Kahler potential, which will cor- 
rect the Kahler potential at higher orders in the inverse volume expansion. As discussed 
above, these are not derived from string theory but are included in a phenomenological 
fashion to ensure an inflection point in the volume direction giving inflation at small vol- 
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umes. In this respect the detailed form of these terms is not important - any other terms 
producing similar physics are equally acceptable. 

We first explain the relation between the model of ( |3.1[ ) and the 1-modulus example 
discussed in the previous section. We shall work at large volume, in the regime of parameter 
space that applies during the runaway period, and start by dropping the C and D terms of 
equations ( |3.1| ), which rapidly become irrelevant at large values of the volume. We also drop 
terms suppressed by higher orders in volume. It is in this region that the 1-modulus model 
is a precise limiting case of the 2-modulus potential. The supergravity scalar potential for 
the large-volume model of (|3.l| ) is 

_ 8^- s (a s A s ) 2 e- 2a ^ 4W a s A s r s e- a ^ £W 2 6 

3V v 2 2V 3 V 2 1 ' 

We have included an uplift term (parametrized by 5) to ensure the physical minimum has 
vanishing cosmological constant. 

There are two fields, t s and t&. These two fields have very different characteristic mass 
scales. Working in the vicinity of e~ asTs ~ V -1 , we have (with Mp = 1) [2] 

m 2 ~K s ~ s —~ — m 2 ~K bl —~— (3 3) 

s h 

As the scale of the potential is V ~ H 2 ~ pj, it follows that m Ts 3> m Tb ~ H. We stress 
that these relations hold in the large-volume regime of the potential, which will correspond 
to the runaway period after inflation has ended. In this regime it is therefore possible to 
integrate out the heavy t s field to generate a single-field potential for the volume modulus 
Tft. This approximation will work with increasing accuracy as the volume increases, as the 
parametric separation between the masses of r s and r& becomes increasingly large. Now, 

<9t, 2V V Va.oTo 



2^JV V V ln V 

We can use equation ( |3.4|) to integrate out t s . Up to terms subleading in In V, this generates 
a potential 

4W 2 (lnV) 3 / 2 jW 2 5 
3a 3/2 V 3 + 2V3 + V 2 ' 

Relating V ~ ~^/2 T b ^° ^ usm S ^ = V 3/2 log(r;,) we see that this potential indeed takes 
the form of the 1-modulus potential ( |2.1| ). The C and D terms generate higher corrections 
of order y~ 1( V 3 and V~ n / 3 , as appropriate to match onto ( |2.1| ). We emphasise that at 
smaller values of the volume, and in particular during inflation, the above procedure of 
integration out is not valid. In this regime the two potentials are qualitatively similar - in 
both cases there is an inflection point in the volume direction - but it is not the case that 
the one-modulus model arises as a strict limit of the 2-modulus model. 
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Having explained the relation of the 2-modulus model to the simpler 1-modulus case 
studied in the previous section, we now turn to the full two-field system described by ([O]). 
We use the following numerical parameters 

£ = 9, C = -173.405, D = 1200, W = -0.1, A = 1, a s = —. 

The choice of numerical parameters is such that a minimum exists at exponentially large 
volumes with a TeV gravitino mass, while the energy scale during inflation is V ~ 10~ 17 Mp. 
As this model has various phenomenological difficulties, to be explained further below, we 
shall not attempt a detailed comparison with observation. The parameters C and D 
are chosen such that at smaller values of the volume there exists an inflection point in 
the volume direction. The large flux numbers present in flux compactifications (typically 
f G3 A G3 ~ x(M) ~ 500) make the relatively large values of C and D not unreasonable. 
We choose the initial conditions so that fields start near the inflection point, where infla- 
tion will occur, ending with the volume rolling away from the inflection point towards a 
decompactification limit. 



lua 




Figure 7: The evolution of the scale factor with time. 



In addition to the volume direction, there is also the r, direction. As the corrections 



to the Kahler potential of (3.1) depend only on the volume, they do not induce a potential 
for the t s direction at constant volume. This direction is lifted only by non-perturbative 
superpotential terms and so is flat at large t s . This is equivalent to the method used 
to generate a flat inflationary direction in Kahler moduli inflation [9,10]. In the vicinity 
of an inflection point for the volume direction, and with the t s direction lifted only by 
non-perturbative effects, the resulting potential is therefore flat and suitable for slow-roll 
inflation. 5 

We choose initial conditions in which the fields have no initial velocities and have initial 
values 

T bMU = 4440, T s , init = 23.56962555 
Figure shows the evolution of the scale factor and figure || shows the evolution of the 

5 There may be additional loop corrections in r s that would lift this flatness, but these are beyond the 
scope of this paper. 
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Figure 8: The evolution of the tj, and t s fields with time. The end of inflation is marked by 
runaway in the t\, direction and oscillations in the r s direction. 

fields with time. For numerical reasons when performing the numerical integration we have 
restricted to the last ten efolds, but there is no conceptual difficulty in extending this to 



field trajectory, followed by a period of oscillations. The initial conditions used are chosen 
to ensure that as inflation ends there is runaway in the volume direction (the direction) 
and oscillations in the r s direction. The presence of oscillations ensures that some quantity 
of post-inflationary radiation will be generated by the decays of the oscillating field. The 
large fine-tuning of initial conditions is primarily due to the need to ensure oscillations, 
and thus radiation, are present at the end of inflation, as it is the radiation that allows 
the attractor solution to exist and be located. In the appendix we show that the tracker 
solution present in the 1-modulus case is also present in the 2-modulus case. Numerically, 
we can confirm that this is also an attractor. This radiation can act as a seed allowing the 
attractor solution to be located. It is the tracker solution that will guide the fields to the 
global minimum. 

In figures |8] and ^ we see that with a purely classical evolution the magnitude of 
oscillations of r s grows with time. If the magnitude of oscillations is not sufficiently reduced 
by particle decays, then in fact the oscillatory r s field escapes from its stabilized location. 
In this case there is no way to find the global minimum of the scalar potential and the 
solution will decompactify. Particle decays are intrinsically quantum events and are not 
explicitly included in the numerical evolution shown in figures ^. If we assume a brane 
wrapping the cycle T s , then the t s direction couples to radiation through the Lagrangian 
terms 



sixty efolds. In figure || we show the form of the field evolution. Inflation occurs early in the 



C = XtsF^F^ + K s - s d^ s d^T s + V(t 8 ) 



(3.5) 



From this the decay width for r s — > 77 can be computed [47] 




(3.6) 
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Figure 9: A plot of the field evolution, showing the inflationary era followed by oscillations in the 
r s direction. 



Evaluating this for the model of 3.1, we have 



- (In V)M P M P 
K ~V, m Ts , — 

The most important relation is the relative size of the decay width T and the Hubble scale 
H. The volume scaling of this is given by 

IV^ 77 _ A 2 (lnV) 3 
H VV 

The question of whether the decays reduce the oscillation amplitude sufficiently rapidly 
is determined by the model-dependent details of the coupling of r s to radiation and on 
the initial conditions. If necessary it is always possible to tune the initial conditions to 
ensure there are sufficient oscillations of the t s field in order to generate sufficient decays 
to radiation. 

The most important point that this model illustrates is that inflation ends with run- 
away in the volume direction together with a small quantity of radiation. The radiation 
will serve as a seed for the attractor solution that will guide the fields into the global min- 
imum. In the specific case of the model above the radiation is provided by decays of t s 
oscillations. 



4. Conclusions and Challenges 

Let us discuss what has been done in this paper and what remains to be done. This paper 
has proposed a mechanism in which high-scale inflation and low scale supersymmetry can 
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coexist. The model discussed explicitly is by no means problem- free and should be regarded 
primarily as an illustration of the main features of the proposed mechanism. These involve 
a runaway at the end of inflation, with the true vacuum far away from that applicable 
during inflation. 

Let us discuss some open questions for this explicit model. First, volume modulus 
inflation was obtained by fine tuning the parameters C and D to generate an inflection 
point. The amount of fine tuning is at a similar order to that required in models of brane- 
antibrane inflation but is not explicit as it relies on higher order quantum corrections to the 
Kahler potential that are not controlled. The use of the volume modulus for the inflaton 
was only for simplicity and it is an interesting challenge to improve the implementation of 
inflation. One interesting generalization would be to a hybrid model of inflation in which 
the decompactification direction that generates the runaway plays the role of the tachyonic 
field that ends inflation. 

A second issue relates to the generation of radiation and isocurvature perturbations. 
In the inflationary regime, we had the inflaton (volume) modulus light to drive inflation but 
also required the second (small) modulus to be light. We also needed a large fine-tuning 
of the initial conditions. This was to ensure that the small modulus would oscillate at the 
end of inflation and generate through its decays the radiation necessary for the attractor 
solution. Indeed, if the second field is very heavy, it always follows the time-dependent 
position of its potential determined by the slowly-changing inflaton. In this case it will not 
experience oscillations after the end of inflationary stage, and so it would not be possible 
to generate the radiation required for the attractor solution. 

The lightness of this modulus is a source of isocurvature perturbations that are highly 
constrained observationally. For the purposes of this paper, which aims simply at combining 
high scale inflation with TeV super symmetry, this is less of an issue. To make our model 
fully realistic, one may consider the second field to be a curvaton and later on having the 
isocurvature perturbations converted to adiabatic perturbations [40]. Another possibility 
to use the second light field for generation of adiabatic perturbations is discussed in [41]. 
In both cases, inflation may produce perturbations of metric with a measurable degree of 
non-gaussianity. This can be a potential advantage of our model, having in mind the recent 
controversy with respect to the possible non-gaussianity in the WMAP data [43-46]. 

On the other hand, the fact that both moduli are light in the inflation regime could 
be modified by additional quantum corrections to the Kahler potential. Such corrections 
are potentially a problem for Kahler moduli inflation as proposed in [9]. In that case the 
flatness of the potential for the small modulus is natural up to and including the leading 
a' corrections to the Kahler potential. String loop corrections [42] not included in [9] have 
the potential to lift that flatness. A correction of the form 5K ~ l/(Vr s ) would give a 
dominant contribution to the mass of r s , resurrecting the r] problem for Kahler moduli 
inflation but avoiding the second light direction in our case. 

All the above issues refer specifically to the particular inflationary model described in 
this paper. In general, for the proposed mechanism to occur the principal requirements on 
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the inflationary model are only that at the end of inflation there is both a small amount of 
radiation present and runaway in the volume direction. In this respect the overall scenario 
is not tied to the particular model used in this paper. The requirement of radiation at 
the end of inflation is suggestive of models of brane inflation, where the brane annihilation 
that ends inflation will generate large quantities of radiation. 

Finally, as generally holds for models of cosmology with TeV-scale supersymmetry, 
this scenario still suffers from the cosmological moduli problem [15]. This was recently 
re-analyzed in [47] for the large volume class of models. While the small moduli decay very 
rapidly, the volume modulus has an MeV-scale mass and gravitational strength interactions. 
It naturally decays very late in the history of the universe ruining successful nucleosynthesis. 
An often discussed solution of this problem is to have a late period of inflation to dilute 
the modulus, such as thermal inflation [48]. Even though an explicit realization of thermal 
inflation in our model is not at hand, the attractor solution automatically contains a high 
fraction of radiation at the end of the runaway period, which may provide the right initial 
conditions for thermal inflation. 

Thus we see that in principle there exist various possibilities to avoid the bound H < 
m 3/ / 2 i n string cosmology. One of these possibilities is to study inflation in KL models [3], or 
in the model proposed recently in [26] . Yet another possibility is to construct a cosmological 
scenario along the lines outlined in this paper. However, none of these possibilities is easy 
and all of them require fine-tuning. This may just mean that the constraint H < m 3 / 2 is 
quite robust and hard to avoid, so one should consider inflationary models where inflation 
occurs on a very small energy scale. 

It is worth emphasizing the fact that low-energy supersymmetry, especially within 
string constructions, tends to put very strong constraints on viable cosmological scenarios. 
The constraint H < m 3 / 2 is just one of them. Other constraints include the cosmological 
problems associated with light gravitino and moduli fields [13,14]. These cosmological 
difficulties may lead us to reconsider how much weight to give to the supersymmetric 
solution of the hierarchy problem. As mentioned in the Introduction, one route out is 
to have models of high scale supersymmetry breaking. These automatically avoid the 
cosmological problems associated to the moduli fields and the gravitino but replace the 
standard solution to the hierarchy problem by appeals to anthropic considerations and 
string theory landscape. On the other hand we can take a positive view on the existence 
of so many cosmological constraints: they provide powerful guidelines for the properties of 
realistic models and eliminate many potential vacua that would be otherwise consistent. 
Fortunately, experimental evidence will soon tell us if low-energy supersymmetry is realized 
in nature in a long anticipated way. 

If supersymmetry is indeed discovered at the LHC, it would falsify many proposed 
models of inflation in both string and field theory. It would then be of extreme impor- 
tance to understand how low-scale supersymmetry can be consistent with the observational 
evidence for inflation at very high energy scales Vi n f » 10 11 GeV, particularly if tensor per- 
turbations were observed by either WMAP or Planck. In this paper we have proposed that 
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this can occur if inflation ends with runaway towards a decompactification direction. The 
energy scales during inflation and the energy scales in vacuo can then be hierarchically 
different. As long as some radiation is present or generated at the end of inflation, an 
attractor solution can guide the fields into the global minimum of the potential, avoiding 
the cosmological overshoot problem. We have presented a model to illustrate this scenario, 
where the volume modulus was the inflaton and the runaway direction was towards the 
large volume minimum. The ultimate aim of model-building is to construct a model that 
is fully realistic when confronted with both cosmological and phenomenological data - our 
work is one step towards this goal. 
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A: Tracker solution for two modulus model 

In section |2| we studied the tracker solution in the presence of radiation for the single field 
potential that arises from integrating out the heavy r s mode. Here we show the existence 
of a tracker solution in the presence of radiation for the full two-modulus model with non- 
canonical kinetic terms. The tracker solution will apply after the end of inflation, namely 
during the runaway regime of the potential. The model is described by 

V = -^-j= (rl /2 - t s 3/2 ) , T b = T b + ic b , T s = t s + ic s . 

K = -2hx(V + £), W = W + Ae- asT % 

V = e K (k^ DiWDjW - 3\W\ 2 ^ . 

Here we have used simply the potential arising for the large volume models, without the 
additional terms C and D that were only important at small volumes. Starting with a 
Lagrangian 

C = j d A x^R + K i3 d^d^ + V{$, (4.1) 

the equations of motion for the complex fields are 

+ 3H& + T) k &$ k + K fj djV = 0, (4.2) 
H 2 -\ (/Qj<W + V + p 7 ) = 0, (4.3) 
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with p 7 the background (radiation) energy density and r§ = IC nl ^f. We fix the axionic 
parts of T b and T s and write the equations solely in terms of the real parts and t s . The 
equations of motion are then 

t* + ZHt 1 + T) k f j i k + -K fj d T , V = 0. (4.4) 

2 

It is useful to rewrite the Priedmann equation ( ^4,3| ) as 

ij = -Kfjfifj - ~p T (4.5) 



In a tracker solution the kinetic, potential and background energies are all constant fractions 
of the total energy. This implies p 1 oc H 2 , K^fifj oc H 2 , and V(t) oc -fT 2 , giving 

H oc -F 2 , (4.6) 

and so in any tracker solution 

The assumption of a radiation background, for which p 7 ~ a~ 4 , then implies a ~ 1 2 and 
iT = i. As the magnitude of the potential is determined in terms of the compactification 
volume by V ~ V~ 3 ~ r 6 9//2 , we can also deduce that 

V ~t 2 / 3 , r 6 = ^ 4 / 9 . 

As at its minimum the t s direction has characteristic mass m Ta 3> H, we anticipate that in a 
tracker solution d Ts V ~ <C V. We will show below that this assumption is self-consistent. 
From equation ( |3.4[ ) this implies 

3 2 

r s ~ In V ~ — In 17, = — In t + constant + . . . (4-8) 
2 3 

We can now check that the above ansatz is consistent and indeed leads to a solution of the 
equations of motion. 

Considering the leading terms at large volume, the Kahler potential K = — 21n(r^ //2 — 
t s ) gives the following Kahler metric 



3 


-9v^ 








3 




a / — 3 / 2 



2 



K=\ JZr. *T h ^=1 3 8^ I- (4-9) 



s. 

4r s r b 



From this we can compute the Christoffel symbols 

1 bb — ! 1 bs — ~ < 1 



V s — ^ Ts r s — r s ~^ r s ~ ^ 
4r 2> ife - i ^- 4 r b ' ss ~4r s " 
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Using the above ansatze we evaluate 



2 

n + 3fff 6 + fljtiTj = -At~ u/9 + (subleading in t), (4.10) 



1 (At s 1 1 



t s + 3Hf s + FfjiiTj = ^ ^— + - - —j + (subleading in t), (4.11) 

We expect any tracker solution to be valid during the regime in which the potential is 
dominated by the third term (the a' 3 correction). During this regime, d Tb V ~ I* 1 the 
tracker solution, we also expect to have d Ts V ~ (as the heavy t s field should be fixed at 
its minimum). In this case, 

l -K b ~id T] V = l -(K bl d Tb V + K bs d Ts V) 



We have here used the assumption that d Ta V <C V. Comparison with equations (| 
and ( [4. 10| ) then gives 



27 t 2 



The t s equations of motion give 



±K s ~id Tj V = ±(K s ~ b d Tb V + K sS d Ts V) 



4 I — 3/2 

-9r 5 U + VT f b (d Ts U) (4.13) 



Comparison with equation (f£4) and (4.11) then gives 



A A I — 3 / 2 

|^ = 9r s y-^pM^) (4.14) 



Using V = 2^5-, eq. ( 4.14 ) is satisfied so long as 

d ^ = ^ L ^r- 4.15 

Equation ( [4.15 ) can be solved to determine the precise value of t s in the tracker solution. 

As this corresponds at leading order in t to our ansatz (|4.8|), this shows that our 
approximations were self-consistent. As the characteristic scale of d Ts V is V itself, the 
requirement of ( 4.15| ) that d Ts V ~ K <C V is consistent with the interpretation that t s has 
been integrated out. 

It is interesting to analyze the origin of the kinetic energy in the tracker solution. 
Writing 

Kijfifj = Kbbfbh + 2Kb s fbf s + K ss f s f s 

= 27^ + 0(t " 8/3)+ ° (t ~ 8/3) ' (416) 
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we see that the kinetic energy of the solution is dominated by the motion of the light field 
(t&) and that the contributions of the heavy field t s to the kinetic energy vanish at large 
t. We can verify that the above results are consistent with the single-field attractor. We 
have 

2 4 1 

V= — x, KE= — 5 , H=—, 
Tit 2 Tit 2 ' Tt ' 

and so = ^ , Q,k e = §f , ^7 = • This reproduces the attractor values of eq (f2.S| ) 
for the single-field evolution. It is thus possible to explicitly solve the equations of motion 
for the 2-modulus model for the runaway regime in the presence of radiation, finding as 
expected the same tracker solution as seen when the heavy modulus is integrated out. 

The tracker solution requires the presence of radiation - in the absence of radiation, the 
small field cannot be confined within its trough and the moduli evolution enters a kination 
phase which eventually leads to decompactification. Furthermore, the simple existence 
of the tracker solution does not imply that the moduli evolution will locate the tracker 
solution. To locate the tracker solution, sufficient radiation needs to be generated at the 
end of inflation to attract the moduli evolution into the tracker solution. Whether sufficient 
radiation can be generated or not is a dynamical question that depends on the amplitude 
of oscillations of the small field and on its couplings to matter - we saw in figure I that m 
a purely classical evolution, with no quantum decays, the small field cannot be confined 
within its trough. 

As for the 1-modulus case, the tracker solution is (numerically) also found to be an 
attractor solution. If the fields can locate the attractor solution, then this will guide them 
through the runaway epoch and into the large-volume minimum without overshooting. 
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